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On ' Abstract 



Analytic formulae are presented for the two-loop perturbative QCD corrections to 6 ^ c 
decay at the zero recoil point, which are required for the extraction of |Vbc| from measure- 
Mh. ments of exclusive B -^ D*lv decays. The results are in agreement with those in g, ^. Some 

comments on the numerical evaluation of the diagrams involved are made. 



1 Introduction 



X 

5t , The magnitude of the Kobayashi-Maskawa matrix element Vcb can be determined experimentally 

by observing the decays of B mesons produced in e'^e~ collisions. In particular, one method 
requires the measurement of the rate of the exclusive semileptonic decay B — > D*lv (see, eg., 
[Q]). From the decay rate, measured at the zero recoil point, i.e., in the special kinematical 
situation where the D* is produced at rest in the B rest frame, one can then extract the value of 
\Vcb\- This method has been used by experiments where B mesons are produced on the T(4S') 
resonance Q and on the Z resonance Q. The statistical and systematic errors are currently of 
the order of 5%. However, in the future, with the CESR collider at Cornell running at increased 
luminosity, the asymmetric B factories at KEK and SLAG coming into operation, and further 
S-physics experiments to be conducted at the hadron accelerators Tevatron, HERA and LHG 
(see, e.g., jQ]), the errors are expected to come down to around 1% or less [g]. Thus, it is 
important that the theoretical input, that is needed to extract \Vcb\ from the measured decay 
rates, be known with equal precision. The purpose of this paper is to present one component 
of that theoretical input, namely the second order perturbative QCD corrections to the decay 
of a 6 quark into a c quark at zero recoil. The other part of the theoretical input consists of 
non-perturbative corrections which are described by an expansion in the heavy quark masses 
rrib and rric- For a review of heavy quark theory and further references, see, e.g., H. 



The reason for using the zero recoil point is that at that point, the non-perturbative contri- 
butions are suppressed by a factor of Ag^^/m^, because of an additional symmetry that exists 
in the infinite quark mass limit. At the same time, however, the zero recoil condition simplifies 
the kinematics of the decay b -^ cW to such an extent, that a complete analytical calculation 
of the perturbative corrections at the two-loop level becomes feasible. 

At tree level, the amplitude for b — > cW is proportional to u{c)'j^ (1 — 75) u{b). In higher 
orders, this gets modified into u{c)j^ (r/v — 7/^75)^(6), where tjv^a are given by perturbation 
series in as- (Only ija enters the expression for B — > D*lv.) At zero recoil, no other Lorentz 
structures appear. The order Ug contributions were calculated in 0. The Feynman diagrams 
that contribute to 'qy^A in order a1 are shown in figure ||. They were calculated analytically 
in ^], confirming the results of a Taylor series expansion in [mi, — m^jmy, [^ that had been 
obtained earlier by one of the authors of Q. In this paper, we present an independent analytic 
calculation of r\y^A- Although our results are expressed in a slightly different way, they are 
completely equivalent to the results of ^. Thus, we confirm the conclusions of ||8|, ^. 

Section ^ describes the main steps of the calculation. Although the details are different, the 
methods we have used are nevertheless related to the ones employed and extensively discussed 
in Q. The major part of the work is the calculation of a set of scalar integrals, many of which 
contain infrared divergences. One example is looked at more closely, and various ways we checked 
our calculation are discussed. The final results are presented in section ^ and are followed by 
our conclusions in section H As a by-product of our work, we developed a new version of a 
numerical technique for evaluating a class of two- loop Feynman diagrams [^] . The new version 
is much more suitable for dealing with some of the rather special diagrams that occur in this 
calculation. However, because the method itself is quite general and may be useful for other 
problems as well, it is explained briefiy in an appendix. 

2 Calculation of rw^A 

An important consequence of the zero recoil condition is that there is no phase space available 
for gluon bremsstrahlung, b — > cW g{g)^ and thus, only the virtual corrections shown in figure |^ 
are needed. We calculate them in the Feynman gauge and use dimensional regularization for 
both ultraviolet and infrared divergences. We neglect the masses of all the light quarks, and we 
leave out the diagram where a t quark loop is inserted into the gluon propagator. With these 
restrictions, all diagrams can be written in terms of the following nine propagator denominators: 



Pi = (/ + A;).(2pi + / + A;) , P3 = /■(2pi + /) , P5 = ^.{2^1 + fc) , Pi = k^. P9 = {I + k) 
P2 = {l + k).{2p2 + l + k), P^ = l.{2p2 + l), Pe = k.{2p2 + k) , Ps = I 



(1) 



where k and / are loop momenta and pi and p2 are the four-momenta of the incoming b quark 
and the outgoing c quark, respectively. Normally, seven of the denominators Pi would be 
linearly independent, but because of the zero recoil condition, which means that pi and p2 are 
proportional to each other, m2Pi = rnip2, we have two additional linear relationships between 
them: 

miPi - 7712^3 = (mi - m2)Ps, ttiiPq - 7712^5 = {mi - m2)P7- (2) 

This greatly simplifies the calculation. 

After projecting the diagrams onto the form factors 7/y,yiQ and expressing all scalar products 
of k, I, pi and p2 that appear in the numerators in terms of the Pi, which we do using Form |jll|], 

^ In the case of r]A, we use an anticommuting 75. 
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Figure 1: Irreducible Feynman diagrams contributing to 6 ^ cW at order a^. The dotted line 
in diagram e represents a Faddeev-Popov ghost. The fermion in the loop in diagram / can be 
either a light quark, a 6 or a c quark. 



we obtain for each one a combination of scalar integrals of the form: 

d'k d'l p.,pj p., , (3) 

where d = 4 — 2e (some of the Ui can be negative). The scalar integrals that occm' can be 
classified into three groups. The first group consists of integrals that can be factorized as a 
product of two one-loop integrals, e.g.: 



// 



d*^^ ({^i — I (^^i X [ d^k . (4) 



The second group consists of integrals in which either z^2 = z^4 = z^6 = or z^i = f 3 = z^5 = 0, 
so that they only depend on one mass scale. Such integrals also occur in on-shell fermion self- 
energies and in anomalous magnetic moments, and can be calculated using recurrence relations 
based on integration by parts [12, 13]. A detailed explanation of the algorithm is given in |14]. 



The third group contains the non-factorizable graphs that depend on mi and on m2- Al- 
though one can still derive relations between them using integration by parts, those relations 
are more complicated than in the case of just one mass scale. We find them to be quite useful 
in a few cases, but are still left with rather a large number of cases we have to calculate from 
scratch. 

Before actually calculating the remaining integrals in the third group, it is worth while to 
investigate their analytic properties by solving the Landau equations. One finds that they can 
have singularities when mi = or ?Ti2 = 0. At mi = m2 there are no singularities because, 
although some poles in the propagators coincide at that point, the integration contours are not 
pinched. However, there may be singularities at mi = m2 on the analytic continuation to higher 
Riemann surfaces. The following very simple one-loop example illustrates these properties: 

If we analytically continue the right hand side of this equation in, say, m2, going around the 
branch point at m2 = and then back to m2 = m-i, log(?7i2) changes into log(m2) + n(27ri), 
which no longer cancels log(?TT.i), and as a result, a pole appears at m-i — m2 = 0. In addition to 
the singularities just mentioned, a certain subclass of the two-loop integrals (^) can also have 
singularities when they are analytically continued to the point mi +m2 = 0. The integrals in this 
subclass correspond to graphs that can be cut into two pieces by removing exactly three massive 
quark propagators. They only occur in diagrams ci, C2, and C3 in figure ||, and in diagram /, 
when the quark in the loop is massive. While this discussion of analytic continuations to negative 
masses and higher Riemann surfaces may seem academic, keeping these properties in mind while 
actually doing the integrations can be very helpful, because it tells us which polylogarithms we 
can expect to occur in the answer. 

We have done most of the integrals needed using Feynman parametric representations of the 
kind described in [^]. Sometimes, it is convenient to differentiate a diagram with respect to one 
(or both) of the masses first, in order to reduce the number of different kinds of propagators, and 
thus the number of Feynman parameters, and then reintegrate with respect to the mass to get 
the final answer. For the integration constant, we can take the equal mass point, which belongs 
to the second group discussed above. Another reason why we might want to differentiate with 
respect to masses, is to make an integral less infrared divergent. 

Let us take the six-propagator integral 








Figure 2: The scalar integrals Ii . . . Iq. The momentum pi enters from the left, p2 —Pi enters at 
the vertex marked 0, p2 leaves at the right. The thin (thick) lines symbolize quark propagators 
with mass mi (777-2) • The dotted lines denote massless propagators. A line with a heavy dot on 
it means the corresponding propagator is squared. 



shown in figure |2| as an example. It is one of the contributions that come from diagram C2 in 
figure |l[ This particular integral has no ultraviolet divergence, but it does contain an infrared 
divergence coming from the region where both k and / are small. Our goal is to express it in terms 
of other integrals in which one (or both) of the massless propagators P7 and Pg are cancelled. 
Let pi^2 = fni,2<l where q is a. fixed four-vector with q^ = 1. First, consider differentiation with 
respect to m2. Using the identities. 
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{{mi - 7772 ) 7772 /l(7n,l, 7772)} = ("7l - 7772) /2("T'l , "^2) " ■"1-2-^3 ("T-l, "^2) , 



12(7711,7772) 



73(7771,7772) 



d'^k d'^l 



d'^k d"/ 



1 



PlP2P3PiPs ' 
1 



(7) 
(8) 

(9) 

(10) 
fll) 



PiPsPePjPs 

The integral I2 is already finite, but it can be simplified further by repeating the above procedure 
with 7771 instead of 7772: 

d 



dm,] 



{(7772 - 77^1)7711/2(7771,7712)} = (7772 - 7711)14(7771,7712) - 7771/5(7711,7772) 



himi,m2) 
75(7711,7712) 



d'^k d^/ 



d'^A; d^/ 



1 



PiP2PiPi ' 

1 
P^PsPiPs ■ 



(12) 

(13) 
(14) 



On the other hand, I3 is stih just as divergent as Ii. However, since it no longer contains Pi, it 
is now easy to cancel Pg: 






^— {mi/3(mi,m2)} =/6(mi,m2) = J J d°'A: d"'/ p2p2p p ■ (1^) 



From figure |2|, it is obvious that I^{mi,m2) = I^{m2,'rni). Thus, in order to obtain the integral 
Ii, all one needs to know are the two finite, four-propagator integrals I/^ and Is, which are 
relatively easy to do by Feynman parametrization, and, as an integration constant, the value of 
13(772,771), which can be calculated by recurrence relations. Neglecting terms of order e, we find 

13(777,7^) = vr^r2(l + e)!!^^|i + 3C(2)-4|, (16) 

14(7771,7^2) = -^|3C(2)--(Li2(u)-Li2(-tx)) + ^^^logfi±^)|, (17) 

where the dimensionless variable u is defined as (7711 — m2)/{mi + 7772), and, after performing 
the final integrations over ttt-i and 7112, 

11(7771, 7772) = n'r\l+e) '^fp,'^ \- + 3C(2) - 12 + - (u2 (n) - Li2 (-n) + log (\^))] ■ 

4:171^1712 '-^ ^V \1 — ^//J 

(19) 
It is comforting to see that the result is symmetric in 7771 and 70,2, as it clearly should be. 

We will not present formulae for all the non-factorizable two-mass scalar integrals (|3|) that 
appear in the diagrams of figure |l], but mention that they can all be expressed in terms of the 
following basic set of logarithms and polylogarithms: 

log(i±^), log(l + «), Li2(n), Li2(-n), U2 {^) , ^^^^ 

Lis (^) , Li3 (^j , Li3 (^j , Li3 (^j , Li3 ((^j^) , Lis [j0T)^) ■ 

These functions are all real and analytic in the range — 1 < u < 1, which corresponds to real, 
positive masses r77i, 7772. It is easy to verify that their only singularities are located at u = ±1, 
and in some cases u = 00, in accordance with the general properties inferred above from the 
Landau equations. The functions that have branch points aX, u = 00 only occur in diagrams 
with massive three-particle cuts, and the trilogarithms only occur in five-propagator integrals. 
We applied various checks on the results for the scalar integrals. First of all, we evaluated 
the ones that are finite in d = 4 by direct numerical integration in momentum space, basically by 
a straightforward extension of the method originally proposed in [|l^] for the two-point two-loop 
"master" diagram with general masses. However, there are a few finite integrals that cannot 
easily be evaluated by this method, unless an important modification is made, which is explained 
in the appendix. Then, there are a number of "sunset" integrals, which we checked numerically 



by means of dispersion relations [16|, and a few other two-point functions we compared with the 
program package Xloops p7| ]. 

A simple analytical check that can be applied to all two-mass integrals is to compare the 
first few terms of their Taylor expansions in u with what one gets by expanding the propagators 
around the equal mass point (7771 = 7772) under the integral sign (as was done from the outset 
in [^). While this is a powerful test, it is not sensitive to possible mistakes in the integration 

6 



constants in diagrams calculated by differentiating and then reintegrating with respect to the 
masses. Therefore, for some integrals, we also looked at their asymptotic behaviour in the limit 
when mi or m2 vanishes (u — > =Fl)- In this limit, the diagrams are again reduced to one-scale 
integrals, but it is more complicated than the equal mass limit because, in general, the diagrams 
become more strongly divergent when one of the masses vanishes. Nevertheless, it is possible 
to obtain the correct asymptotic behaviour by expanding under the integral sign, provided one 
adds certain counterterms corresponding to the Taylor expansion of certain subdiagrams of the 
diagram considered, following the prescription for asymptotic expansions on the mass shell given 
mH]. 

To conclude this list of consistency checks, we mention two tests that we applied to the 
contributions of complete diagrams, rather than individual scalar integrals. We verified that if 
mi and m,2 are interchanged, diagrams ai and 02 in figure |l| are transformed into each other, and 
similarly 61 <-> 63, 62 <-*■ &2j ci <-> C3, C2 <-> C2, (i <-> d, e <-> e and / <-^ /. This is a consequence 
of charge conjugation symmetry. Secondly, by putting the masses of the two external quarks 
equal to each other, we reproduced the two-loop on-shell quark wave-function renormalization 
constants in [^] . Except for one part — the contribution^ of a quark loop whose mass is neither 
zero, nor equal to that of the external quark line — which we calculated separately, we could 
obtain all the wave-function renormalization factors by simply substituting u = in our formulae 
for contributions to r]v- 

3 Results 

Here, we present our results for r]v,A up to 0{a^) in QCD with N colours and Nl flavours of 
massless quarks, but no t quarks. The colour factors are given by 

iV^ - 1 1 

Ca = N, Cf = , Tp = - . (21) 

Following fs], |9[, we renormalize the strong coupling constant in the MS scheme at the scale 
M = ^mhrric but use on-shell renormalization for the masses. This renormalization procedure 
has the nice property that it respects the symmetry under m^ <-^ rric mentioned at the end of 
the previous section. Therefore, i]v,a are even functions of the variable: 

TJih — rrir 

u = — . 22 

mb + rric 

This symmetry was exploited in [Q by using the variable p = {nib — rric)'^ / (mbmc) = 4ti^/(l — -u^) 
rather than 6 = {rrib — mc)/rnb = 2n/(l + u) as an expansion parameter, in order to obtain a 
more rapidly converging series. 

It turns out that the trilogarithms Lia ( -^^ j and Lis ( :j^ ) 1 which appear in the contribu- 
tions of diagrams ci, C2 and C3, cancel out in the sum ci -|- C2 -|- C3. The other functions ( pO|) 
appear in the following combinations: 

(23) 
(24) 
(25) 



i = 


io.(;+") 

\l-uj 


Ci = 


Li2 (u) — Li2 {—u) 


C2 = 


Li (" 'M 1 'e' 


^''[u + lj ' 4^ 


^Equation (22) in [O 


1 






Li3 ^-7 + L13 + ^^^ + 1:^ L12 —— 

,u + l/ \u—\} 6 3 V^i+1 



Lis 



Au 



(n+l)2 



+ Li3 



-4ii 



\u 



1)2 



+ ^^^ + ?^Li. 



2u 
ii + 1 



-C(2) (41og(l + tz) - 2^) - -£ (Li2 (u) - Li2 [-u]) . 



(26) 



(27) 



Note that i, Ci and £2 are odd functions of u, while £3 and £4 are even. The abbreviations a 
and Z3 are defined by 



as(v/'^fe'^c) 
C' = ^ 

zs = C(3)-41og(2)C(2) 



(28) 
(29) 



Finally, the results are 
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1 + aCp {£ --6 



+a^ 



CfTfNl 



u 3 



(-104n2 + 464^"^ + 1144^^ + 32u^ 
(u2 - 1)4 

,2 (-16 + 224n2 + 896^"^ + 416n6 + 16ii^) 
(n2 - 1)4 
32 + 128ii2 + 1616u^ + 1312^6 + 48u^) 



+CfTf { C(2) 



+ . 



.^ + 200^2+344^4) 



u{v? - If 



u{v? - lY 



,6 



128 I 296„,2 I 344„,4' 



+c, 



+Cl C(2) i 



32 - 128u^ - 2368U'' - 3200u'' - 480u» -^ + ^u^ + ^u 



u(n2 - 1)^ 



iu' - 1) 



(-16) 



+ C(2) 



(-32 + 48^2) i_m 



J-16^ 

+^2^(^2_l)+^3 



n2-l 
24 - 96^2 53 ' 

-n + — 



+ 



+e 



u 



+CFiCA - 2Cf) C(2) i 



-8u) 



n2-l 



+ C(2) 



10^2 



17 

u 



- 1^2 + n^4 

ii2(ii2 _ ]^) 



2 + 2n2 

^2-1 



4u (-16n) 12-6^2 (_g^2) 17- 

+^-1^5 7 + i-2^ ; ^41,3 + La^-12 ^ 1- ^3—2 



u^-l 
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u 



u^-l 



u^-l 



u^-l 



1 +aCF 

W CfTfNl 

+Cpr^|c(2) - 



10^ 
3 ' 



u 



+ T 



-64 - f ^2 + 2|2^4 + 2|2^6 + M^8) 
(u2 - 1)4 



32, 



^ ^ (-f + 120n2 + '-fu') ^ ^2 f + ¥^^^ + ^^^ + ¥^" + f ^° 



m(i(2 - 1)2 



(w2 - 1)4 



(30) 



^'^^ n(u2 - 1)4 

32 _ 1664 2 _ 8000,, 4 _ 7808,, 6 _ 992 8 1016 , 584 2 , 704 4 ' 
3 3"- 3» 3" 3" 9 "TgU'T'g"' 

+-'-'2 r"2 TV4 '" 



16 _ 32^2 (_8-64n2) 



+/:2 1 o ':^ +J^: 



2 




' u(n2 — 1) u 

( 32 _ 56„,2 40 _ 46„,2 , 70„,4 61 

+CACA - 20 c(2) i ^^^ + c(2) ' ,f ;:,^^ + ^ ^ 

(-1 + f ^2 + 1^4 + 4^6) I + 1^2 + ^^4 



V?{v? — 1)2 u(li2 — 1)2 

(31) 



.A '-^-^'''-^-''' -24.3.Ai±J^..3^-'°"' ^"^ 



m(u2 - 1)2 ^ ' ^ ^,2 _ 1 ' •^ ^2 _ 1 9 

Note that every single term in the expressions (|3^) and (^) is manifestly analytic in u for at 
least \u\ < 1. This fact, combined with the symmetry u <-> — n, proves (once more) that rjv^A 
can be represented by convergent power series in u^ for all finite positive values of mi, and rric- 

4 Conclusion 

We have performed an independent calculation of the vector and axial vector form factors rjv 
and rjA that describe the decay of a 6 quark into a c quark at zero recoil up to order a'^. We have 
compared our expressions ( |30| ) and (31) with the corresponding formulae in [^ by rewriting the 



latter in terms of the functions (|20D used in this paper, and found that the two calculations are 
in perfect agreement. Both confirm the series expansion in the mass difference mb — rric obtained 
earlier in [Q. 

The actual numerical values of rrih and rric are such that a few terms of the series expansion 
are already sufficient to achieve the accuracy that is needed in practice {v? ~ 0.29). On the 
other hand, the exact analytical formulae are not only valid for b ^ c decays, but also for other 
cases, such as t ^ 6, where the mass ratio is larger, and the series expansion converges more 
slowly. They also allow one to study the limit when one of the quark masses goes to zero, which 
is impossible if only a limited number of terms of an expansion in the mass difference are known. 

In the appendix, we have identified the cause of a numerical difficulty that arises when the 
two-dimensional numerical integration method of [|^ is applied to certain two-loop diagrams 
with several coinciding thresholds, and suggested an alternative method that solves the problem. 

We would like to thank J.G. Korner for suggesting this project and P. Post for helping us 
check some sunset diagrams. We are also grateful to A. Czarnecki, A.G. Grozin, K. Melnikov, 
M. Neubert and D. Pirjol for discussions. This work was partly supported by the Graduiertenkol- 
leg "Element arteilchenphysik bei mittleren und hohen Energien" at the Johannes Gutenberg 
University in Mainz. 

Appendix 

In this appendix, some details of the numerical integration methods we used to check the scalar 
two-loop integrals appearing in this work are given. We restrict ourselves to cases that are both 



ultraviolet and infrared finite in d = 4. Basically, we follow the approach proposed in ||T^ for 
the numerical integration of the scalar two-point two-loop master diagram with general masses. 
The fact that the diagrams occurring here are three-point diagrams is not a problem: since the 
external momenta pi and p2 are proportional to each other, they have a common rest frame, 
and all the steps in the derivation of the integral representation given in |jl^ can be repeated 
with only a few trivial changes in the formulae. However, the fact that pi and p2 are on their 
respective mass shells does give rise to some difficulties, which we discuss here. 




\ a 



-2m^ -2m2 



Figure 3: Singularities in the complex /cQ-plane originating from the quark propagators P5 and 
Pq (C3), and from the gluon propagator Pj (0). While the integration contour can always be 
rotated away from the real axis by an angle a, it can only be shifted away from the origin in 
diagrams not containing the gluon propagator Pj. 



By writing the loop momenta k and / in the common rest frame of pi and p2 as {k(),k±), 
{loJ±), and integrating out their space components k±, l±, we obtain representations of the 
form 

/■OO /"OO 

dkodloF{ko,lo;p) (32) 



for the scalar integrals (|3|). Here, we have explicitly indicated the dependence on the small 
imaginary quantity ip which is added to all propagator denominators in accordance with the 
causal prescription. As a function of the complex variable ko, the integrand F{kQ,lo;p) has 
singularities (branch points) located at the positions shown in figure |3[ Identical pictures can be 
drawn for the complex variable Iq and for the combination feo + ^o- When the limit p — > is taken, 
the singularities move onto the real k^ and Iq axes. In general, this makes the integrals along the 
real axes more difficult to evaluate numerically and in some cases even divergent. Sometimes, 
for example when the scalar integral under consideration does not contain either of the gluon 
propagators Pj and Pg, the problem can be solved by rotating the integration contours away from 
the real axes by an angle a (which must be the same for ko and Iq), Iq -^ e*'^/o; ^0 ~^ e*°A;o, and 
subsequently shifting (one of) them away from the origin, e^^ko — > e*"fco~/3, as shown in figure ^. 
Scalar integrals we checked in this way include I2 in figure |2| and J J d kd l/^PiP^P^Ps). 

This escape is not possible in diagrams where both the ko and Iq contours are trapped at 
the origin by gluon propagators, such as, e.g., 



d^k dH 



1 



PiP-iPePrPs 



(33) 



We stress that if the limit p — > is taken after all integrations have been performed, one obtains 
a finite result for the integral (|33|)P|. However, one would like to set p to zero, or at least, 
very close to zero, before performing the /cq and ^0 integrations in d32), and, if one attempts 



In the case of equal masses mi = m2, it reduces to the one difficuh one-scale integral A'^(l, 1, 1, 1, 1) 
6^(2) log(2) — I (^(3) needed in two-loop QCD or QED corrections to on-shell fermion propagators |1£|]. 
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to do so, a non-integrable singularity appears at the point k^ = Iq = 0. More specifically, 
F{Xko,Xlo', P = 0) ^ 1/A^ as A ^ 0. In the original integral (p3[), this logarithmic divergence 
conies from the region where k± and l± tend to zero, while /cq and /q are of order k"^, tj_, as one 
can see by rescaling 



k± -^ Xk± , l± -^ Xl± , k(j — > A^^o , ^0 -^ ^'^^0 ■ (34) 

Under this transformation, the integration measure in ( |3^ ) scales like A^'', while, for small A, 
the integrand goes like 1/X^^. 

The problem can be solved very easily by interchanging the order of the integrations. That 
is, we first perform the fco, Iq and all angular integrations analytically, leaving two integrations 
over k± = \k_i_\ and l± = \l±\ to be done numerically. 

Below, we shall give an explicit formula for the resulting integrand. Although our main 
reason for deriving it was the need for an independent check on our analytic results for (|33| ) and 
a few other, equally nasty cases, the formula is valid for a much more general diagram: 

where 

Di = {k+ pi)^ -ml + ip Di = {1+ p^)"^ -m\ + ip 

D2 = {k + p2f -ml + ip D^ = {l+p^,f -ml + ip (36) 

D3 = {k + l + p-sf -ml + ip. 

with the restriction that all momenta pi are proportional to one another. Apart from that, the 
Pi and the masses rrii are arbitrary. Depending on the choice of the pi, (^5|) then corresponds to 
a diagram with two, three or four external legs. We use this notation for the sake of flexibility, 
even though it has some redundancy (which could be used, for example, to set P2 = P5 = 0). 
Working in the rest frame of the momenta pi, we find the following representation: 

J = 4vrM dfeW d^M^^^|^ ^^"'^"^-^j + "^-^°^ + ^3-(-^i-P?-^4 




with 



^^12^^45 ^^12^^45 



+ (1 ^ 2) + (4 ^ 5) + (1 ^ 2, 4 ^ 5) L (37) 

FH.) = >o,(^40). (38) 

C^ = {±u,-p^+p^f-u'j, (39) 



X + p|] ± ng 


; 


p^+p^f- 


2 


U4 


= 


M5 


= 



and 

ui = \Jm\ + k'j^ — ip U4 = yml + l^ — ip 

U2 = yml + k\-ip Ms = \Jml + l\-ip (40) 

tif = yrnl + (k^ it l^Y — ip . 
The representation (|37|) can be applied directly as it stands to the case (|33|). In general though, 



LliC C"^ 

avoided by a rotation of the contours: l± — > e~*"/^, fcj 



there will be poles near the real axes at the points where the Cj • vanish. They can easily be 
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